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Nonequilibrium Pump-Probe Photoexcitation 
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Relaxation of electrons in a Hubbard model coupled to a dissipative bosonic bath is studied to simulate the pump- 
probe photoemission measurement. From this insight, we propose an experimental method of eliciting unoccupied part 
of the single-particle spectra at the equilibrium of doped-Mott insulators. We reveal first that effective temperatures of 
distribution functions and electronic spectra are different during the relaxation, which makes the frequently employed 
thermalization picture inappropriate. Contrary to the conventional analysis, we show that the unoccupied spectra at 
equilibrium can be detected as the states that relax faster. 


1. Introduction 

Time-resolved and out-of-equilibrium dynamics of elec¬ 
trons in pico to femtosecond time scale has been re¬ 
vealed by rapidly developing pump-probe photoemission 
spectroscopy.^ A challenge is to get insight into emerging 
properties of correlated electrons from the relaxing electrons 
far from equilibrium. 

Instead of equilibrium Green’s functions used to analyze 
the single-particle excitations from the equilibrium,^ nonequi¬ 
librium Green’s functions have been used with a finite reso¬ 
lution limited by laser pulse width.^ However, understanding 
the full relaxation process including long-time scale remains 
a challenge. 

Electron relaxations were often analysed phenomenologi¬ 
cally as if they were in equilibrium at each time steps but at 
time-dependent temperatures different from the state before 
pumping^"^ such as two-temperature model'^. However, we 
will show later that these quasithermal treatments are insuf¬ 
ficient. 

In addition, the quasithermal approaches do not tell us 
microscopic information such as so-called dark sides (un¬ 
occupied sides) of the spectra above the Fermi level.In 
fact, clarifying the momentum-resolved dark-side spectra is 
claimed to hold a key for understanding challenging issues 
such as the mechanism of the pseudogap and high-Tc super¬ 
conductivity in the cuprates.So far, the dark side has few 
accessibility such as the inverse photoemission with severely 
limited energy resolution. 

In this paper, by introducing a simulation scheme for 
nonequilibrium dynamics of many-body electrons, we first 
show that effective temperatures determined for the electronic 
spectra, and the distribution (occupation) of electrons follow 
time evolution clearly distinct each other, which unjustifies 
the quasithermal model. In addition, we in fact find that elec¬ 
tron relaxations strongly depend on excitations beyond the qu- 
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asithermal picture. This inhomogeneous relaxation emerges 
in a way that simpler electron-hole excitations relax faster. We 
propose an experimental method by which unoccupied (dark 
side) quasiparticle states in equilibrium are identified by faster 
relaxation after pumping. 

2. Model and Method 

2. 1 Hamiltonian 

Here, we simulate relaxation processes of photo-excited 
states of doped Mott insulatorsby including phase relax¬ 
ation process in addition to energy, and momentum relax¬ 
ations.^^ We introduce a minimal model that consists of a fi¬ 
nite Hubbard ring described by Hw coupled to bosons, and 
2-level atoms to understand the generic feature of strongly 
correlated electrons in nonequilibrium. 

The energy and momentum relaxations of electrons in 
crystalline solids occur through electron-phonon interactions. 
Simulating phonons carrying finite momenta are necessary 
for the relaxation. The phase relaxation is more complicated: 
In a closed system, the phase relaxation of the system never 
occurs. To simulate the phase relaxation, we introduce addi¬ 
tional two level atoms as detailed below. At short intervals, 
we couple them to and disentangle them from the system con¬ 
sisting of electrons and phonons. When we disentangle these 
atoms, we introduce random phases and mimic the phase re¬ 
laxation. 

Our hamiltonian reads 

H{t) = ^h( 0 + + ^fb + ^jc, (1) 

which is schematically illustrated in Fig. 1. A Hubbard ring 
under a time-dependent electric field dFfQjdt is defined as 

HuV) = - j F + h.c.] + %, (2) 

i,o- i 

where, W is the band width, (c.^) is the operator creating 
(destroying) an electron with cr-spin at the i-th site, and U 
is the on-site Coulomb repulsion. Here the electron number 
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Fig. 1. (color online): Schematic description of our hamiltonian. A Hub¬ 
bard ring Hwif) coupled to phonons a+q^ is excited by a time-dependent elec¬ 
tric field Fit). These phonons are also coupled to two-level atoms. We ad¬ 
ditionally introduce phonon dissipation with a phenomenological decay con¬ 
stant 7G and with randomizing the phases of wave function (p (see Appendix 
A.l and A.2 for details). 


operators are defined as h.^ = ^Icr^io-' 

We introduce bosons that exchange energy and momentum 
to realize the energy, and momentum relaxations, which are 
subtantiated by (acoustic) phonons in real solids. Two modes 
of bosons (denoted by v and jj) are introduced by the Hamil¬ 
tonian H\y described by 

Hb = HoJb Yj ’ ( 3 ) 

where al (aj is the creation (annihilation) operator of a boson 
mode with the label u = y,ju with energy Here v = qo 
and Id = -qo specify momenta ±qo = ±27x1 L, where L is the 
number of sites. The coupling of the bosons to the electrons 
are described by i^fb! 

^ I] («/. + + h-c.], (4) 

k=27T£IL 

where g is the coupling constant. These bosons change the 
total momentum of the electron wave function with minimal 
increment and decrement in momenta for the L-site Hubbard 
ring. We also introduce a 2-level atom coupled to each mode 
of the boson. These atoms absorb and supply the bosons. The 
coupling between the two 2-level atoms (labeled by a and JS) 
and bosons is described by as 

^ u=v,M 


photo-excitation and relaxation processes of the correlated 
electron systems will be studied in the following part of the 
present paper. Before going into details, we remark the gener¬ 
ality of the present one-dimensional model as a prototype of 
strongly correlated electron systems including those in higher 
dimensional systems. One may suspect that the Tomonaga- 
Luttinger liquid behaviors may govern dynamical processes 
of the one-dimensional Hubbard ring which should be in a 
strict sense distinct from the universal Landau’s Fermi liquid 
behaviors in higher dimensional correlated metals. However, 
finite-size systems with finite time span and finite energy res¬ 
olution in the present study do not distinguish the Tomonaga- 
Luttinger liquid behaviors from those of the Fermi liquid and 
do not hamper us from extracting general tendency of the dy¬ 
namical processes in the correlated electron systems. 

The relaxation of the dark-side spectra is proposed to be 
dominated by the relaxation of simple particle-hole excita¬ 
tions. The particle-hole excitations in the one-dimensional 
electron systems hardly remain as well-defined qausiparti- 
cles under the presence of the electron-electron interactions: 
The particle-hole excitations in one-dimensional interacting 
electron systems, for example, described by the Hubbard 
ring are reconstructed into collective excitations orthogonal 
to fermionic excitations in the thermodynamic limit, which 
is characterized by the Tomonaga-Luttinger liquid behaviors. 
However, the quasiparticle residue Z remains finite even in 
the one-dimensional systems provided that the Hubbard ring 
with finite sites are examined: The quasiparticle residue Z in 
the one-dimensional interacting electron systems is scaled by 
the power of the system size although Z approaches zero as 
the system size increases.Therefore, even though the one¬ 
dimensional Hubbard ring is employed, the relaxation pro¬ 
cesses in the numerical simulation given in the following sec¬ 
tions are not specific to the one-dimensional systems except 
the relatively slow spin relaxations nearly decoupled from the 
charge excitations. 

2.2 Method for time evolution 

The wave function of the Hubbard ring coupled to the 
bath can be expanded by using the basis representing the N- 
electron eigenstates of the Hubbard ring |0^^), the Fock states 
of the phonons \ny) and |^^), where Uy and are numbers of 
phonons in the each modes, and the eigenstates of the two- 
level atoms \a) and \P), as given as, 

Adim Tb-1 Lb-1 

l<i>(0> = EZEEZ Cm\ny,nn\afiif) 

m=l ny=0 ni^=0 a=a,bp=a,b 


where is a diagonal Pauli matrix and are lad¬ 
der operators interchanging the eigenstates of hK is the 

coupling constant between the 2-level atoms and bosons. Here 
the energy difference of the levels between the ground state 
and the excited state of the atom, denoted by \b} and \a), re¬ 
spectively, is chosen as the energy of the boson mode hco^. 

By employing the time-dependent hamiltonian H(t), the 


x|<D^)(8iK)0|n^)®|Q'>(8i^>. (6) 


As detailed below, the time-evolution based on the wave 
function is twofold. Together with the ordinary time-ordered 
hamiltonian dynamics 


|O(0) = T exp 



dtH{t) 


mto)}, 


(7) 


2 
















J. Phys. Soc. Jpn. 


FULL PAPERS 


Started with the ground state wavefunction |0(ro)), energy, 
momentum and phase relaxations of the total system are phe¬ 
nomenologically achieved through the stochastic surrogate 
hamiltonian approach^^ and a damping constant yc of the 
boson numbers,as detailed in Appendix A and B. The 
former realizes repeating and successive decoupling of the 2- 
level atoms from the bosons, and reset of the 2-level atoms 
to their ground states accompanied by randomization of their 
phases, everytime after a fixed interval of the time evolution 
by the Hamiltonan H. 

We use a specific form of F(t), 

F{t) = f cos [copuiseC?' - to)], (8) 

JtQ 

to represent applied pulsed lasers with the amplitude Aq, the 
pulse center to, the pulse width td, and the frequency u^puise- 

2.3 Spectral functions 

We calculate the momentum k and frequency oj dependent 
spectral weight A{k,(jj\ t) and its occupied part Aocdk, oj’, t) 
from the time evolution, by projecting the bosons and the 
two-level atoms out and by taking a long-time approximation 
that we call stretched-time Lehmann-representation methods. 
The stretched-time Lehmann representation is nothing but a 
Lehmann representation based on a projected electronic wave 
function |O(0)proj as detailed in Appendix C. The projected 
wave function is prepared by operating a projection operator 

^0 = {fj 

m 

(g) {ny= 0| (8) {a = b\® {j3 = b\, (9) 

to the wave function |O(0) as 

Aditn 

|O(0>p,oj = ^ol®(0> = L (10) 

m=l 

where a^iit) = Cm-ofi-^bit). 

This method enables the frequency representation 
of transient phenomena although it satisfies Freericks- 
Krishnamurthy-Pruschke’s resolution limit only approxi¬ 
mately because of the finite probe-pulse width.^ Typically 
0(10^) fs pulses give us the spectrum resolution ~ 0(10“^) 
eV, which justifies our results given below. 

2.4 Effective temperatures 

Here we show how to estimate the effective temperatures 
of the physical quantities. Our estimation is based on results 
of the angle-integrated spectra, 

AoCci.^’> 0 — AQQQ^k, (jJ',t), (11) 

k 

and 

A{oj,t) = YjA{k,oj;t). (12) 

k 


The distribution function is defined as 

n{oj, t) = Aocc(^, t)IA{(jj, t), (13) 

which corresponds to the Fermi-Dirac distribution if an effec¬ 
tive temperature becomes well-defined. 

When we fit the time-dependent spectra by an exact elec¬ 
tronic spectra at a finite temperature Tqq, we can estimate ef¬ 
fective temperatures from two different functions, namely, the 
temperature of nico, t) and that of A^oj, t). Here, we introduce 
the exact electronic (occupied) spectra at finite temperatures 
A^®^(6 g), keTeff) and the distribution functions 
Then we can estimate the temperatures estimated from the 
following function by minimizing 

/](?; t,«) = f'doj ^eff) - f(co, i)f. 

JAi 

(14) 

The effective temperature that minimizes 



(15) 


(16) 


are denoted by and TQfij^{t), respectively. If the photo- 

excited states are really described by the thermally excited 
ones, these two temperatures should be the same, namely, 
TtWnif) = TeffA(0- 

3. Results 

We simulate relaxation of photo-excited strongly correlated 
metals. The spectra of the strongly correlated metals in equi¬ 
librium show a typical structure regardless of dimensionality 
of the systems: The spectra consist of the incoherent upper 
Hubbard band located around oj - U, the incoherent lower 
Hubbard band located for oj < 0, and the low-energy coher¬ 
ent band around the Fermi energy While the upper 

Hubbard continuum is well separated from the low-energy 
coherent band, the latter and the lower Hubbard band touch 
each other or are not separable, in general. For the later dis¬ 
cussion, we call the unoccupied part of the low-energy exci¬ 
tations including the quasiparticle band as mid-gap states. We 
especially highlight relaxations of dark side mid-gap states of 
doped Mott insulators after photo excitation by a laser pulse. 

Below, we set the energy unit W/4 = fi = 1. Parameters 
used for real-time evolutions are = Ni = 3, L = U = 
8, ojo = V 2 / 4 , g = 0.05, A = tan“^ 197r/40 ~ 0.08, and 
yoAt = 0.05 with a time interval At = 0.1. We use pulses 
with frequencies ct^puise = V2/4 much lower than direct gaps 
between occupied states and upper Hubbard bands. Typically, 
laser pulses used in this paper give energy ~ 0.01 mJ/cm^ to 
electrons if we set W/4 = 0.5 eV and the inter-site distance 
as 3 A. The unit time of the present system corresponds to ~ 
1.32 fs. 

Here we use pulses with frequencies u^puise much lower 
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than direct gaps between occupied states and upper Hubbard 
bands. Even though the photon energy is not enough for direct 
transitions to the upper Hubbard band, both of the mid-gap 
low-energy states just above the Fermi level and the upper 
Hubbard bands are excited due to the short width of pulses 
and the energy-time uncertainty relation. 

3.1 Relaxation of mid-gap states, upper Hubbard bands, 
and effective temperatures 

First, we show relaxation of the density of states in Fig.2. A 
pulsed laser with duration td = 3 and amplitude Aq = 4 excites 
the system and induces nonequilibrium mid-gap excited states 
within an energy window, yUF-0.85<6G)< F’uhb = 4, and 
upper Hubbard bands (F’uhb = 4 < 6l)), where yup is the Fermi 
energy. Here, F’uhb is chosen as an upper bound of the low- 
energy spectrum around the Fermi energy yup. 

We observe that the upper Hubbard band relaxes first. The 
mid-gap states absorb the relaxed excitations from the upper 
Hubbard band and it joins in the relaxation. 



0 500 1000 , 1500 

time t 


Fig. 2. (color online): Relaxation of occupied density of states for f//r = 8 
with 6 electrons {N^ = Nf) on a 8-site ring. A laser pulse with td = 5,Ao = 4, 
and ujpuise = V2/4 is used. The center of the pulse is located at r = 0. Inset 
shows time-dependence of integrated spectral weights for mid-gap states and 
upper Hubbard bands, /mid and /up. A hnite line width <5 = 0.2 is used for the 
spectra. 

By using the following integrated spectral weights, 

^£2 ^+00 

lEuE2if) - I dioAocci^^t)! I d(jjAocci(^^t), (17) 

kJEi 0—cx> 

we here define the integrated spectral weights for the mid¬ 
gap states and upper Hubbard bands as /mid(0 = W^e^ubO) 
and ^up(0 — ^Fuhb ■+oo(t), respectively. The integrated spectral 
weights /mid and are summarized in the inset of Fig. 2. The 
relaxation rate of /mid up to f ~ 10^ depends on the square 
of g when other parameters are fixed, while the relaxation at 
longer time scales is additionally affected by phase relaxation. 
We employ g to reproduce the experimental relaxation rate of 
doped Mott insulators^’ in the time range up to r ~ 10^, ir¬ 
respective of ultrashort-time relaxation mechanisms that may 
depend on dimensionality and details of the hamiltonian.^^"^^ 

The numerical results of the effective temperatures T^f^fft) 
and TeffACO are shown in Fig.3, where the energy windows 



Fig. 3. (color online): Time-dependent effective temperatures of distribu¬ 
tion functions reff„(0 and spectral functions are shown by white cir¬ 

cles with;^“^ and;^“^ in color plots. Inset shows double peak structures of 
Xn^(r, Teff) around t = 740, which induce abrupt changes in 


are chosen to extract low-energy contribution as Ai = -3 and 
A 2 = +3 (Ai = -00 and A 2 = ^uhb) are set for T^f^nit) 
(TeffA(O) in Eq.(14). Around the occupied band bottom and 
the upper Hubbard band distant from the Fermi energy yup, the 
density of states in non-equilibrium tends to deviate from the 
counterpart in equilibrium more significantly than the density 
of states around the Fermi energy yUp (for example, see Fig.l 
ofRef.28). 

The effective temperatures TQf^fft) and TQf{^(t) are signifi¬ 
cantly different as shown in Fig.3, which clearly demonstrates 
that the spectral weights relaxes faster than the distribution 
functions. The deformation of the spectral weights is nothing 
but a manifestation of the electron correlations, which inval¬ 
idates simple effective electron temperature approaches. The 
deformation of the spectral weights also means that the ex¬ 
cited carriers above the Fermi level are not injected to the 
dark-side of the equilibrium spectra, in contrast to a naive ex¬ 
pectation for visualizing the dark-side with hot carriers. 

In contrast to the relaxation of the spectral function, the re¬ 
laxation of the distribution function is not smooth as shown 
in the time dependence of (see Fig.3). The abrupt 

change in the time dependence of is a manifestation 

of a strong deviation from the equilibrium high-temperature 
distribution for a short time after the irradiation, t < 800. 
Similar strong deviations from the equilibrium distribution in 
initial stages of hot carrier relaxation are reported for phonon- 
assisted relaxation of hot carriers in free electron systems.^^ 

3.2 Angle-resolved and unoccupied spectra through loga¬ 
rithmic differences in relaxation processes 

Next, we examine momentum-resolved nonequilibrium 
spectra. In Fig. 4(a), we show the equilibrium momentum- 
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resolved spectrum with the same parameter set used in Fig. 2 
and 3. Just after the photo-excitation, the occupied spectrum 
^occ(k, w; t) becomes completely different from the spectrum 
in equilibrium shown in Fig. 4(a). As evident in Fig. 4(b), the 
peaks in nonequilibrium spectra are not coincident with the 
peaks in equilibrium unoccupied states. 

Therefore, to observe or estimate the dark-side equilibrium 
spectra of doped Mott insulators, photo-excited electrons are 
seemingly useless. However, after closer inspection, we real¬ 
ize that, from the photo-excited spectra, the dark-side equilib¬ 
rium spectra can be estimated as proposed below. 



Fig. 4. (color online): (a) Spectrum of a doped Mott insulator out of 
equilibrium, Aocdk, aj\t) in comparison with total spectrum in equilibrium 
A(k, oj; -oo). The shaded area shows occupied spectrum Aocdk, uj; t) and the 
solid curves show total spectrum A(k, co; -oo) in equilibrium. The shaded area 
appears at energy ranges where no equilibrium unoccupied states exist, (b) 
Time-dependent spectra at smallest unoccupied momentum k = njl from 
t = 10 to r = 950. The nonequilibrium spectra have a prominent peak at the 
equilibrium quasiparticle peak oj = EqpCk = 7:12) (downward arrow) , while 
they have peaks at oj different from the energy of the equilibrium peak as 
well, as indicated by the upward vertical arrow. A finite line width S = 0.2 is 
used as in Fig.2. 

Here, we propose an approach to extract the dark-side spec¬ 
tra in equilibrium, which may replace the standard effective 
temperature approach. We focus on energy- and momentum- 
resolved relaxation rates of Aocdk, oj', t) and frame a hypothe¬ 
sis: The relaxation rates of excited nonequilibrium states have 
local maxima when the states are located at the unoccupied 
quasiparticle energy in equilibrium. The hypothesis seem¬ 
ingly contradicts Landau’s Fermi liquid theory,in which the 
quasiparticles’ life times are longer than incoherent states in 
equilibrium. However, outside the equilibrium, the quasipar¬ 
ticle excitations have efficient relaxation pathways as detailed 
below. 

Then, we explain how the photo-excited states in the dark 


side consisting the quasiparticles relax faster than the inco¬ 
herent excitations. First of all, we assume correspondence of 
a many-body excited state with small total momentum |^ol 
to unoccupied quasiparticle peak due to added electrons in 
equilibrium, which is a central notion of the Fermi liquid the¬ 
ory. In equilibrium, the detailed balance between emission 
and absorption of the phonons and the Pauli blocking guar¬ 
antee long life times of quasiparticles even in the presence 
of the electron-phonon couplings.However, out of equi¬ 
librium, where the detailed balance and the Pauli blocking 
become inefficient, single particle-hole excitation easily de¬ 
cays by emitting energy and momentum through the electron- 
phonon couplings. Due to the short-time violation of the en¬ 
ergy conservation, the particle-hole excitations decay through 
first order perturbation process of the electron-phonon cou¬ 
pling g. Therefore, by recalling the first assumption, out-of- 
equilibrium states injected at the quasiparticle peaks decay 
faster than more complicated excited states such as multiple 
particle-hole excitations. 

To extract spectral components with large relaxation rates, 
we propose and use logarithmic differences of nonequilibrium 
occupied spectra Aocdk, oj; t) for some time interval from t\ 
until t 2 , which gives an averaged relaxation time Ti-eiax(^, ^) 
as 

-1 _ 1 ApcciK co\ ti) 

^relax(^9 Oj) ^2 “ H Apcdk, Oj;ti) 

By this average, the effects of unavoidable noise in experi¬ 
ments are suppressed. The slow relaxation of /mid in compar- 



Fig. 5. (color online): (a) Time-dependences of mid-gap weights /mid- The 
partial swap introduces noise in /mid- Inset shows parameter sets of td and Aq. 
The circle A corresponds to the parameter set td = 5 and Aq = 4 used in Fig. 
2 and 3. (b) Logarithmic differences and Treiax in nonequilibrium occupied 
spectra at the smallest unoccupied momentum k = njl. When laser pulses 
with Aq > 3 are used, the logarithmic differences show significant minima. 
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ison with the probe laser pulse width (-100 fs) validates our 
stretched-time Lehmann representation, where t 2 - h (> 0) 
longer than the relaxation time of /mid is also better. 

In Fig. 5, we show the usefulness of the logarithmic dif¬ 
ferences of the nonequilibrium unoccupied spectra. For suffi¬ 
ciently large amplitude of laser pulses, we obtain significant 
minima of the logarithmic differences as functions of oj. The 
minima agree with the quasiparticle peak, as illustrated in Fig. 
5(b). In contrast, for smaller amplitude of laser pulses, the 
minima become less significant and are not coincident with 
the quasiparticle peak as shown in Fig.D-2 of Appendix D. 

It is crucial to use the laser-pulse amplitude Aq in an appro¬ 
priate range. The laser pulse must generate particle-hole exci¬ 
tations efficiently beyond the linear-response regime in a con¬ 
trolled way, which requires Aq to be comparable with that re¬ 
quired for the Bloch oscillation, ABioch- We propose to choose 
Aq in an extended range roughly given by Aq > u^puisse > ^ ltd 
but not Aq » \ ltd. See the detailed condition for the pump 
laser in the Appendix D. 

4. Summary and Discussion 

In this paper, we have simulated relaxation of photo-excited 
doped Mott insulators by employing a Hubbard ring coupled 
with bosonic degrees of freedom. We have shown the ex¬ 
istence of two distinct effective electronic temperatures that 
make the standard effective-temperature approaches unjusti¬ 
fied. We instead propose a method to extract the dark-side 
equilibrium spectra from nonequilibrium pump-probe relax¬ 
ation by using logarithmic difference of the time-dependent 
spectra. The logarithmic difference offers a stable experimen¬ 
tal tool that circumvents inevitable noises in the experiments. 

We have also explored sufficient conditions for detectable 
signals in experiments: Amplitude of the pumping laser is 
required to satisfy Aq > u^puise > I ltd- The laser pulses 
with a sufficiently large amplitude generates photo-excited 
states consisting of combinations of elementary excitations 
and quasiparticle excitations. It enables us to extract the equi¬ 
librium dark-side spectra from the photo-excited nonequilib¬ 
rium spectra. 

The elementary excitations in equilibrium have relatively 
fast relaxation paths in the photo-excited nonequilibrium 
states. The proposed method opens a novel way to access the 
momentum-resolved dark-side spectra of the strongly corre¬ 
lated electron systems. 

The newly proposed scheme to access the dark-side spec¬ 
tra will be helpful for the studies on unoccupied spectra of 
high-Tc superconducting cuprates. It gives us a direct access 
to structure and formation of the so-called pseudogap in the 
cuprate superconductors, which hold a key for understanding 
the mechanism of the high-temperature superconductivities in 
the cuprates. 
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In this Appendix, we detail the simulation method for 
the energy, momentum and phase realized by the sur¬ 
rogated hamiltonian [17]. The stretched-time Lehmann- 
representation methods are also introduced. Constraint on 
pulse laser for efficient observation of quasiparticle excita¬ 
tions is also detailed. 

Appendix A: Energy, momentum, and phase relaxation 

A. 1 Swapping operator 

As expalined in the main article, we introduce an operator, 
namely, swapping operator that stands for the reset of the 2- 
level atoms to its ground states \h) and replacing the phases 
of the electronic and bosonic wave functions with random 
phases with the periodic time interval Al Our swapping oper¬ 
ator is defined as 

Adim r.B - 1 Lr- I I - 

Dswap|ffi(0) — 

m=l ny=0 n^=0 y a,j3=a,b 

xe'V».;»v.»p loW) (g) |«^) 0 \n^) 0 \b) 0 \b), (A-1) 

where (fm;ny,n^ is a random phase. 

We also use a partial swap defined by the following partial 
swapping operator: 

Dpswap = 1 + dOswap? (A'2) 

where 1 is an identity operator. After applying the Oswap and 
^pswap^ we normalize the wave function. 

A. 2 Additional disspation mechanism 

To control the temperatures of the lattice, which is mim¬ 
icked by the two boson modes, in addition to the swapping 
operators, we use a Landau-Lifshitz and Gilbert-type dissipa¬ 
tion operator to keep the number of bosons small enough par¬ 
tially following Gisin’s idea (see Ref. 18) for the time interval 
At: 

Adim r.B- 1 Lb- 1 

DGisinl^(O) = E2EE2 

m=l ny=0 a=a,bp=a,b 

® \ny) ® \n^} ® |a) ® ^), (A-3) 

where yc is a phenomenological decay constant of the num¬ 
ber of the bosons. Here, we impose the condition that the equi¬ 
librium temperature to be zero, because experimental temper¬ 
ature is normally negligible in terms of the excited energy of 
electrons. 

For practical numerical treatments, the additional dissipa¬ 
tion mechanism with the decay constant yc is helpful since it 
keeps the averaged bosonic number small and, therefore, re- 


6 




J. Phys. Soc. Jpn. 


FULL PAPERS 


duces the truncation errors for the bosonic Hilbert space. In 
a real crystalline solids, the number of the phonon modes is 
large enough to keep the density of the bosons. However, due 
to the limit of computational resources, we employ a minimal 
phonon subsystem. Therefore, we need to get rid of excess 
accumulation of the bosons due to the limited number of the 
phonon modes. The decay constant yc is again useful to avoid 
the accumulation of the bosons. 


Appendix B: Time evolution 

Our time evolution consists of 2 steps as detailed below: 
First, we calculate the hamiltonian time-evolution as 


|0(r/+i)) = T exp 







(B-l) 


rH+i ^(f\ 

where Te is a time-ordered operator. Next, we apply 

the swapping operator. 


1^. XX 

^ = 

y I ^pswap ^Gisin I 1)) | 


(B-2) 


We successively repeat the above 2 steps. The swapping time 
defined by 4wap = 4+i - U (= MAt) controls the phase relax¬ 
ation, and an energy relaxation as well. In this study, we set 
M = 200. 


Appendix C: Ansatz for spectrum 

Here we introduce an approximate method to calculate 
time-dependent spectral functions. The spectral functions 
revealed by time-dependent angle-resolved photoemission 
spectra provide us with physical pictures. The present formal¬ 
ism corresponds to the limit of the long-time window of the 
laser pulse used as the probe (see Ref.3). We call the follow¬ 
ing approximate method as stretched-time Lehmann spectral 
representation. 


Here we apply the following ansatz justified in the steady- 
state limit: Starting from |O(0)proj, we define an approximate 
time-dependent Green’s function for the occupied and unoc¬ 
cupied states as, 

Goccik, w; 0 = |am(0p 

n,m 

^ . o , 77(A-1) Z7(A) ’ 

OJ iS Ej^ — Eyyi 

and 

Guno(^9 ^9 0 ~ ^ l^m(0l 

n,m 

a)+ i5- 

where is the eigenvalue that corresponds to |0^^). 

By using the above ansatz for the Green’s functions, we can 
define the stretched-time Lehmann spectral representation as 

A(k, cl>; t) = --Im [Gocc(^, 0 + Guno(^, 01 • (C4) 

71 

Here we note that the occupied spectrum is given by an occu¬ 
pied Green’s function defined as 

Aocc(^, t) = --Im [Gocc(^, oj\ 0] • (C-5) 

71 

Appendix D: Pumping laser dependence 

There is an optimal parameter range of amplitude and du¬ 
ration of laser pulses for efficient observations of the dark- 
side spectra by utilizing the logarithmic differences of the 
time-dependent angle-resolved photoemission intensities. For 
the efficient observation, the laser pulses are required to ef¬ 
ficiently generate particle-hole excitations including the de¬ 
sired momentum p intended to measure the quasiparticle. Be- 


C.l Projection 

First, we need to eliminate the bosonic degrees of freedom 
and project the entire wave function onto the A-electron par¬ 
tial Hilbert space expanded by as in Eq.(lO): 

A(Jitn 

m=\ 

Here, we expand the projected wave function by the N- 
particle eigenstates of H^iit = 0), {|0^^)}. 

C.2 Ansatz 

As an important observable in experiments, the time- 
dependent spectral functions have often been discussed. How¬ 
ever, the energy-resolved spectral functions are definitely ap¬ 
proximate ones whenever the observations are carried out in 
non-steady states (see Ref.3). If the time-dependence of the 
system is slow enough in comparison with the time scale 
set by resolution limits, the approximate time-dependence of 
the energy-resolved spectral functions becomes well-defined. 


^0 



Fig. D l. Constraints on amplitude Aq and duration of laser pulses to ob¬ 
serve quasiparticles with momentum amplitude \p\. The curves and lines are 
determined by Ao/ujpuise = \pl = const., and td = 1/ujpuise- The dashed 
curve qualitatively illustrates the upper bound of Aq, below which the “adia¬ 
batic transitions” accompanying the particle-hole excitations efficiently occur 
by suppressing the Landau-Zener tunneling. The shaded parameter region is 
suitable for generating a pair of particle-hole excitations. 
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low, we clarify what are the optimal amplitude and duration 
of the laser pulses. 

We start with two trivial constraints. The efficient observa¬ 
tions are hampered in two regions. For In, the laser 

pulse contains less than a single cycle of electric-field oscil¬ 
lations, which is called a subcycle laser pulse. Subcycle laser 
pulses are expected to generate net electric currents strongly 
depending on their pulse shape or phase. Thus, the subcycle 
pulses are not suitable for stable observations of the dark side. 
The other trivial constraint comes from pump fiuence. For 
too large pump fiuence or tdA^ » 1, the system gains too 
much energy. Then, the relaxation processes involves multi¬ 
ple particle-hole excitations and becomes complicated, which 
may be observed as sudden increase in relaxation time as laser 
fiuence is increased. 



Fig. D*2. (a) Time-dependences of mid-gap weights /mid- The partial swap 

introduces noise in /mid - Inset shows parameter sets of td and Aq. (b) Logarith¬ 
mic differences and Treiax in nonequilibrium occupied spectra at the smallest 
unoccupied momentum k = nil. When laser pulses with Aq > 3 are used, the 
logarithmic differences show significant minima, (c) Logarithmic differences 
and Treiax at momentum k = n 12 for Aq < 3. The minima are located at o) 
different from Eqp. Inset shows /mid for weakly excited cases with Aq < 3. 


To gain sufficient intensity of the emitted electrons by the 
probe laser, the pump light is required to generate sufficient 
amount of particle-hole excitations. In the present formalism 


such excitations are generated by the acceleration of electrons 
by the time dependent vector potential F(t)lea, where e is the 
elementary charge and a is the lattice constant. The amplitude 
of the vector potential F(t)lea is controlled by a parameter 
Aq. 

Let us clarify how F(t) oc Aq generates the particle-hole 
excitations. These particle-hole excited states orthogonal to 
the ground state in the electronic subsystem Hu are attained 
by many level crossings generated by the vector potential 
F(t) oc Aq and additionally by the hybridization gap indirectly 
induced by the coupling to bosons, which transforms the sim¬ 
ple level crossing to the avoided crossing. 

When the system originally in the ground state is driven 
by adiabatically passing through such an avoided crossing, 
the state becomes containing particle-hole excitations origi¬ 
nally orthogonal to the ground state in the absence of the light 
and bosons. Namely, the particle-hole excitations are gener¬ 
ated by the “adiabatic transitions” in the driven system. When 
the time dependent vector potential F(t) oc Aq is an oscil¬ 
lating pulse with a frequency u^puise and a gaussian envelop 
exp [-(r - determines the number of times 

that the system passes through the avoided crossing, where 
to is the gaussian pulse center in the time domain and td is the 
pulse width. Therefore, to increase the chance of such “adi¬ 
abatic transitions” to particle-hole excited states, we need to 
drive the system for sufficiently long time td at least to satisfy 
^puhQ^d > 1 provided that the above condition tdAo < 1 is 
satisfied. 

An additional requirement for the “adiabatic transitions” to 
occur effectively exists: If Aq becomes high, the probability of 
the non-adiabatic Landau-Zener tunneling at the avoided level 
crossings increases without following the adiabatic route. It 
prevents effective generation of particle-hole pairs. Therefore, 
for higher Aq, td is required to be longer to achieve the same 
density of particle-hole excitations. 

Equivalently, an appropriate density of particle-hole exci¬ 
tations are obtained in the region (1 - FLz(Ao))^^puise ~ 
where the Landau-Zener tunneling probability 0 < Flz < 1 
(the probability per each crossing through the avoided en¬ 
ergy crossing point) depends on Aq and largely increases 
with Aq. Thus the “adiabatic transition” with the probability 
1 - Flz(Ao) decreases with increasing Aq. This means that the 
upper limit of desired Aq increases with td. 

We also have a constraint that Aq has to be larger than an¬ 
other threshold. By noting that the time-dependent phase F(t) 
causes lattice-momentum shifts, for observation of quasipar¬ 
ticles with momentum amplitude \p\, we need the maximum 
of the time-dependent phase F(t) larger than \p\. For a gaus¬ 
sian pulse, the maximum of the phase F(t) is approximately 
given by AoIoj^uIsq- Therefore, when we observe the quasipar¬ 
ticles with momentum amplitude |p|, we need Aq that satisfies 
Ao/^^puise > \p\- In general, Aq required by the above con¬ 
straint is comparable with that required for the Bloch oscilla¬ 
tion, ABloch- 

By taking the above constraints into account, the suitable 
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parameter region of the pulse laser amplitude Aq and duration 
td is determined for the efficient observation of the dark side. 
The constraints and the obtained parameter region are sum¬ 
marized in Fig.D-1. The appropriate region for the parameter 
values is roughly given by Aq > ~ u^puise ^ ^jtd- 

The parameter values studied in the main text, namely, Aq = 
4, td = 5, u^puise = V 2 / 4 , and \p\ = InjL ~ 0.8 satisfies this 
criteria. 

To demonstrate the general validity of the hypothesis on 
the relaxation rates in the logarithmic differences, we choose 
several sets of parameters td and Aq, in addition to Fig.5 of 
the present paper, and show that the hypothesis works if the 
amplitude of the laser pulse is large enough, here Aq > 3, as 
illustrated in Fig.D-2. 
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